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A  theoretical  and  experimental  study  on  the  feasibility  of  determining 
the  size  of  a  single  submicron  particle  by  observing  its  Brownian  motion 
characteristics  has  been  carried  out.  The  method  is  based  on  measurement  of 
the  particle  motion  interf erometrically  using  the  light  scattered  from  a  pair 
of  intersecting  laser  beams.  The  particle  is  assumed  to  be  in  thermal 
equilibrium  with  a  fluid  medium.  Due  to  the  viscous  damping  of  the  fluid,  the 
motion  of  the  particle  exhibits  relaxation  behavior.  The  relaxation  time  may 
be  obtained  from  the  interference  signal.  Knowing  the  temperature  and  the 
viscosity  of  the  fluid,  the  mass,  hence  the  size  of  the  particle  may  be 
determined  from  the  relaxation  time.  Monte  Carlo  simulation  of  the  Brownian 
motion  detector  response  has  shown  that  that  the  relaxation  time  may  be 
obtained  from  the  statistics  of  the  time  between  extrema  of  the  signal.  A 
proof-of-concept  experiment  was  conducted  using  latex  spheres  of  known  sizes 
as  test  particles.  The  signal  to  noise  ratio  of  the  experimental  setup, 
however,  was  not  adequate  for  determining  the  size  of  the  particles. 


1 .  INTRODUCTION 


1.1  Overview 

Measurement  of  particle  size  In  the  submicron  (0.01  -  0.1  ^m)  range  has 
many  Important  scientific  applications.  For  example,  in  the  study  of  soot 
emission  in  combustors  and  burners,  the  soot  formation  and  oxidation  rates  can 
be  inferred  from  the  soot  particle  size  evolution.  The  combustion  of  coal 
particles  and  the  formation  of  fly  ash  can  be  studied  in  a  similar  manner. 
Submicron  particles  are  also  important  in  corrosion  of  boiler  surfaces  and 
heat  exchangers  in  electric  power  plants  or  chemical  processing  plants  and  in 
erosion  of  turbine  blades  in  advanced  combined  cycle  processors.  Other  areas 
where  an  understanding  of  the  role  of  submicron  particles  is  important  include 
the  nucleatlon  and  condensation  processes  leading  to  acid  rain,  and  soot 
formation  and  enrichment  of  polycyclic  aromatic  hydrocarbons  in  internal 
combustion  engines.  Measurement  of  submicron  particles  is  basic  to  an 
understanding  of  all  of  these  processes. 

In  this  report,  we  describe  an  attempt  to  develop  a  technique  for 
accurate,  in-sltu  determination  of  submicron  particle  size  distributions  based 
on  a  Laser  Doppler  measurement  of  the  Brownian  motion  of  the  particles.  The 
objective  was:  to  review  the  scientific  basis  of  the  technique,  to  conduct  a 
proof-of-concept  feasibility  experiment,  to  determine  through  the  experiment 
the  practical  difficulties  of  the  technique,  and  to  identify  the  solutions 
needed  to  overcome  these  difficulties. 

1 .2  Review  of  Particle  Sizing  Techniques 

In  general,  particle  measurement  in  the  submicron  size  range  can  be 
categorized  as  either  probe  sampling  techniques [ 1 . 1-1  .4 ]  or  in  situ  optical 
techniques .[ 1 .5-1 .13]  In  the  sampling  probe  approach,  a  vacuum  system  of 
varying  degree  of  complexity  is  used  to  draw  the  particle  laden  stream  through 
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grids  or  filter  paper.  The  deposited  particles  are  then  analyzed  by  electron 
microscopy [ 1 .2,  1.3]  or  optical  reflectivity  and  transmissivity 
measurements .[ 1 .4 ]  The  difficulty  in  these  direct  sampling  methods  is  that 
the  surface  deposition  of  the  particles  influence  the  particle  size 
distribution  in  an  unknown  way.  The  physical  presence  of  a  sampling  probe  may 
also  interfere  with  the  phenomena  under  study.  This  interference  is 
especially  severe  when  high  spatial  resolution  is  required. 

In  situ  optical  techniques  have  an  advantage  over  probe  sampling  because 
the  flow  is  not  disturbed.  Optical  techniques  may  be  classified  as  a  single 
particle  counting  or  a  volume  averaged  measurement.  In  single  particle 
counting,  the  measurement  system  examines  only  one  particle  at  a  time  and 
records  the  particle  size.  The  particle  size  distribution  is  obtained  by 
tabulating  the  sizes  of  many  particles.  In  the  volume  averaged  method,  the 
measured  signals  are  a  convolution  of  the  optical  properties  of  the  individual 
particles  and  the  particle  distribution.  Therefore,  various  assumptions  have 
to  be  made  to  obtain  the  size  distribution  from  the  measured  signals. 

The  most  well  developed  optical  techniques  are  based  on  an  analysis  of 
the  amplitude  of  the  scattered  light  from  the  particles.  Measurements  with 
various  optical  configurations  have  been  reported. [1 .5-1 .13]  For  example, 
measurements  and  analyses  were  performed  with  multicolor  systems,  with  light 
of  various  polarizations,  and  on  the  angular  distribution  of  the  scattered 
light.  The  interpretation  of  these  measurements  is  usually  based  on  Mie 
theory. [1 .14,  1.15]  A  comprehensive  review  of  these  measurements  is  given  in 
Ref.  [1.16]. 

The  major  disadvantage  of  the  light  scattering  technique  is  that  the 
scattered  amplitude  is  sensitive  to  the  optical  properties  of  the  particles, 
which  are  in  general  not  well  known.  In  general,  researchers  have  interpreted 
the  experimental  results  based  on  Mie  theory  and  a  refractive  index  for  the 
particles  which  is  based  on  bulk  measurements .[ 1 .4 ,  1.17]  The  validity  of 
using  the  refractive  index  of  the  bulk  material  is  open  to  question. 
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Therefore,  there  is  a  need  to  develop  an  optical  particle  size  measurement 
scheme  which  is  independent  of  the  optical  properties  of  the  particle. 


For  completeness,  we  mention  other  miscellaneous  optical  sizing  schemes. 
These  schemes  include  the  opto-acoustic  cell [1.18,  1.19]  which  is  a  volume 
averaged  measurement  and  is  dependent  on  the  particle  light  absorption 
properties,  and  the  holographic  technique, [1 .20]  which  is  difficult  to  apply 
to  particles  of  size  smaller  than  the  wavelength  of  the  illumination. 

In  this  report,  we  have  studied  submicron  particle  sizing  in  a  gas  stream 
by  measuring  the  Doppler  shift  of  the  scattered  light  from  the  particles  in 
Brownian  motion.  Since  the  scheme  is  a  frequency  measurement,  the  result 
should  be  independent  of  the  optical  properties  of  the  particle.  The  measured 
signal  is  a  function  of  the  gas  viscosity,  which  depends  on  the  temperature 
that  can  be  readily  determined  independently,  and  the  particle  mass,  which  can 
be  related  to  the  particle  size. 

This  concept  of  particle  size  determination  is  not  new. [ 1 .21-1 .24 ] 
Nevertheless,  previous  measurements  based  on  this  concept  have  all  been 
interpreted  in  terms  of  spectral  broadening  of  the  scattered  radiation  due  to 
the  Brownian  motion  of  a  large  number  of  particles.  Therefore,  these 
measurements  were  volume  averaged.  Our  approach,  however,  is  to  operate  in  a 
single  particle  counting  mode  in  which  the  individual  particle  size  is 
inferred  from  the  Brownian  motion  of  a  single  particle. 

1 .3  Organization  of  the  Report 

The  remaining  chapters  of  this  report  are  organized  as  follows.  In 
Chapter  2,  the  basic  physics  of  Brownian  motion  is  reviewed.  In  particular, 
we  have  highlighted  the  statistical  nature  of  the  motion,  which  precludes 
the  measurement  of  particle  velocities  by  the  traditional  laser  Doppler 
technique.  In  Chapter  3,  a  method  for  interpreting  the  information 
contained  in  the  Doppler  signal  of  the  scattered  light  and  a  demonstration  of 
the  concept  in  a  Monte  Carlo  computer  simulation  of  the  signal  are 
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discussed.  The  theoretical  basis  for  the  method  is  described  in  Chapter  4. 

In  Chapter  5,  the  signal  to  noise  required  for  adequate  signal  detection  is 
estimated  and  the  experimental  set  up  for  the  proof-of-concept  demonstration 
is  described.  The  analysis  of  the  signal  obtained  in  this  experiment  is  then 
presented.  It  was  concluded  that  the  signal  is  dominated  by  shot  noise  and 
the  Brownian  motion  information  cannot  be  recovered  from  the  signal.  In 
Chapter  6,  we  examine  some  of  the  physical  limitations  on  the  improvement  of 
S/N  by  increasing  the  beam  intensity.  The  summary  and  conclusions  of  this 
project  are  presented  in  Chapter  7. 
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2.  BASIC  PHYSICS  OF  BROWNIAN  MOTION 


In  this  chapter,  the  basic  physics  of  Brownian  notion  is  described.  The 
nature  of  the  notion  is  such  that  within  the  tine  that  the  particle  traverses 
a  fringe  pair  in  a  laser  Doppler  velocimeter  systen,  its  velocity  nay  have 
changed  drastically.  Therefore  such  a  systen  does  not  yield  meaningful 
velocity  measurements.  The  proper  interpretation  of  the  information  contained 
in  the  Doppler  signal  will  be  discussed  in  Chapter  3. 

2 .1  The  Nature  of  Brownian  Motion 

When  a  small  particle  is  immersed  in  a  fluid  medium,  it  exhibits  an 
irregular  motion  as  a  result  of  the  bombardment  of  the  fluid  molecules.  It 
should  be  noted  that  the  momentum  transfer  to  the  particle  in  each  collision 
is  extremely  small,  but  the  collision  frequency  is  extremely  high 
(~1014/s).  The  motion  of  the  particle  is  the  aggregate  effect  of  a  large 
number  of  collisions. 

A  schematic  representation  of  the  time  dependent  motion  of  a  0.1  ^m 
particle  is  shown  in  Figure  2.1  in  order  to  indicate  the  time  scales 
characteristic  of  the  motion.  In  this  figure,  random  changes  in  particle 
velocity,  of  order  10-b  m»s-1,  occur  in  characteristic  times  of  10-14  seconds 
due  to  collisions  with  gas  molecules.  This  time  scale  is  much  faster  than 
can  be  resolved  by  practical  instruments.  Therefore,  only  an  aggregate 
velocity  is  observable.  This  aggregate  velocity  obeys  a  Maxwellian  velocity 
distribution.  The  characteristic  time  for  the  Maxwellian  distribution  to 
develop  is  the  "relaxation  time"  trelax  8  Figure  2.1),  which  is  a 

function  of  particle  size  and  fluid  properties.  Clearly  if  one  wants  to 
obtain  any  information  on  the  Brownian  motion  of  an  individual  particle,  one 
must  sample  the  time  history  of  particle  motion  in  an  interval  short  compared 


to  the  relaxation  time,  or  the  effects  will  be  averaged  out.  Conversely,  one 
cannot  monitor  the  instantaneous  Brownian  velocity  of  the  particle,  because 
electronics  and  signal-to-noise  limitations  constrain  realizable  sampling 
intervals  to  be  much  longer  than  the  10~i4  s  characteristic  collision  time. 
The  measurement  strategy,  then,  is  to  sample  with  an  interval  on  the  order  of 

.  Q 

10  s,  which  is  fast  compared  to  the  relaxation  time.  The  type  of  signal 
measured  with  this  kind  of  time  resolution  will  be  a  superposition  of  many 
random  events,  but  the  time  dependence  of  this  signal  should  exhibit 
statistical  properties  representative  of  the  Brownian  motion  of  a  particle  of 
a  particular  size. 

The  statistical  analysis  of  the  Brownian  motion  of  a  particle  in  thermal 
equilibrium  with  a  fluid  medium  is  described  by  the  Langevln  equation: [2.1 ] 


dii 

—  -  pu  +  A(t) 


(2.1) 


where  ii  is  the  velocity  and  p  is  the  damping  coefficient.  The  random 
excitation  A(t)  has  the  property  that  the  solution  of  Eq.  (1),  which  is 
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has  the  Maxwellian  velocity  distribution  in  the  limit  of  large  t.  The  time 
scale  for  approaching  this  limit  is 


'  relax 
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In  continuum  flow  (X-nfp  «  d),  for  a  spherical  particle  of  diameter  d  and 
mass  m,  p  is  given  by 
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where  p,  is  Che  viscosity  of  Che  fluid  medium.  For  Che  free  flow  limit 
(^mfp  »  d),  the  value  of  f3  is  given  by  the  Epstein  formula, 
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where  c  is  the  mean  thermal  speed  of  the  fluid  molecules  and  a  is  the 
accommodation  coefficient.  A  value  of  a  ■  0.5  was  used  in  the  simulation  to 
be  described  in  Chapter  3.  Values  for  the  relaxation  time  are  shown  in 
Figure  2.2. 
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As  an  example,  we  consider  a  particle  of  0.2  pm  diameter  with  a  specific 
gravity  of  2.  For  air  at  300  K.  and  at  1  atmosphere  pressure,  the  relaxation 
time  is  "Crgiax  *  150  ns.  To  obtain  information  about  the  velocity 
distribution,  the  sampling  time  should  be  less  than  Treiax.  A  sampling 
time  At  of  the  order  of  10  ns  was  used  in  the  simulation  to  be  described  in 
Chapter  3.  Furthermore,  a  time  between  samples  of  10  ns  was  assumed. 

The  time  scale  for  the  fluctuating  excitation  A(t)  is  of  the  order  of  the 
collision  time  tc  between  the  fluid  molecules  and  the  particle,  which  is 
~ 10  s.  During  each  collision,  the  particle  changes  its  velocity  by 


8kTm 

m 


(2.6) 


where  m  is  the  mass  of  a  fluid  molecule.  The  value  for  oV  is  ~10  m*s 


Since  the  sampling  time  At  is  large  compared  to  the  collision  time  xc, 
the  velocity  sampled  is  the  average  over  At  of  the  time  dependent  velocity 
given  by  Eq.  (2.2).  The  requirement  that  the  velocity  distribution  has  the 
asymptotic  limit  of  a  Maxwellian  distribution  requires  the  sampled  velocity 
Wjj  at  nAt  to  have  a  probability  distribution  of,  [2.1] 
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PARTICLE  DENSITY  -  2  x  10*  kg  m“* 
AIR,  300  K 


Figure  2.2.  Characteristic  Relaxation  Time  for 
Brownian  Motion. 
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Therefore,  If  wn  ■  (un,vn,wn),  and  if  Che  Doppler  system  is  sensitive 
Co  the  velocity  component  v,  Che  probability  distribuCion  function  for  wn  is 
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2 .2  Measurement  of  Brownian  Motion  With  a  LDV  System 

Consider  a  cross-beam  laser  Doppler  veloclmeter  system,  operating  in 
fringe  mode, (2.2]  set  up  to  measure  the  velocity  of  a  particle  in  Brownian 
motion.  For  a  perfect  detection  system  as  a  whole,  the  transfer  function  is 
shown  in  Figure  2.3.  The  transfer  function  is  the  relationship  between  the 
input  quantity  of  Interest,  the  position  of  the  particle  in  Brownian  motion, 
and  the  output  signal,  which  is  the  beat  amplitude  of  the  two  fringe  forming 
laser  beams  as  seen  by  a  square- law  detector. 
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signal  “  <jRe(E^e  +  — 2e  (2.9) 

~  1  +  cos  (Ak«jr) 

In  the  normal  LDV  application,  the  particle  traverses  the  illuminated  volume 
with  a  constant  velocity  (from  A  to  B  in  Figure  2.3).  If  u  is  the  velocity 
component  in  the  direction  of  k^ -^ , 
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Che  signal  will  be  modulated  fully  (A'B')  as 
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(2-11) 
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and  che  velocity  may  be  obtained  from  the  period  x  of  the  modulation, 
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u  -  2*/Ak  x  (2.12) 

Note  that  implicit  In  this  velocity  measurement  is  a  length  scale  (Ak)-1  so 
that  a  time  measurement  (x)  can  yield  a  velocity. 

In  the  detection  of  Brownian  motion,  however,  u  is  a  statistical  quantity 
which  changes  rapidly  in  time.  In  particular,  the  particle  may  change 
velocity  appreciably  (or  may  even  reverse  its  direction)  before  it  traverses  a 
"fringe"  of  the  transfer  function.  On  the  average,  the  "excursion"  of  the 
particle  is 


<z> 

rms 


relax 


/(kT/m) 


(2.13) 


The  values  of  <z>rma/\  are  shown  in  Figure  2.4.  The  average  excursion  is 
smaller  than  the  fringe  spacing  2x/Ak,  therefore  the  signal  is  no  longer  fully 
modulated  and  (Ak)-1  ceases  to  be  a  relevant  length  scale.  Because  of  the 
nonlinear  nature  of  the  transfer  function,  the  amplitude  of  the  signal  will 
not  carry  meaningful  information.  For  example,  when  the  particle  oscillates 
about  CD  in  Figure  2.3,  the  signal  response  C'D'  is  quite  different  from  that 
(E'F')  of  a  similar  oscillation  EF. 

The  only  meaningful  information  contained  in  the  signal  is  the  time 
between  the  extrema  of  the  signal  modulation.  Each  extremum  of  the  signal 
represents  a  reversal  in  direction  of  the  particle,  and  therefore,  a  zero  of 
the  Brownian  velocity  u.  In  order  to  interpret  the  signal,  it  is  necessary  to 
examine  the  statistical  distribution  of  the  zeros  of  u.  This  procedure  will 
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PARTICLE  DIAMETER  (^m) 


Mean  Excursion  Distance,  Compared  with  the  Wavelength  of 
Visible  Light,  for  Particles  (Specific  Gravity  **  2)  in 
Brownian  Motion  in  a  Gas  at  T  •  300  K,  p  «  1  Atmosphere. 
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be  described  in  Chapter  3.  From  the  statistics,  the  relaxation  time  of  the 
particle  may  be  determined,  and  then  the  size  information  may  be  determined 
from  Eq.  (2.4)  or  (2.5). 
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3.  SIMULATION  OF  THE  BROWNIAN  MOTION  DOPPLER  SIGNAL  AND  SIGNAL 
PROCESSING  STRATEGY 

In  this  chapter,  a  numerical  simulation  of  the  Doppler  veloclmeter  signal 
from  a  particle  in  Brownian  motion  is  reported.  The  simulation  was  based  on 
a  Monte  Carlo  technique  for  calculating  the  motion  of  the  particle  from  the 
stochastic  solution  of  the  Langevin  equation  [Eq.  (2-8) j.  This  motion  was 
then  coupled  into  the  response  of  a  laser  Doppler  veloclmeter  system  operating 
in  fringe  mode  to  produce  the  Doppler  signal.  Then  the  results  from  using 
different  methods  for  processing  the  signal  were  compared  to  formulate  a 
signal  processing  strategy  for  the  determination  of  the  particle  relaxation 
time. 

3*1  Simulation  of  the  Brownian  Velocity 

The  measurement  approach  utilizes  a  laser  interferometric  system  for 
measurement  of  an  interference  signal  which  is  sensitive  to  the  particle 
displacement.  Such  a  system  is  shown  schematically  in  Figure  3.1.  This 
approach,  in  simple  terms,  is  to  measure  the  time-dependent  interference 
signal  arising  from  a  single  submicron  particle  (£  0.1  ^m  diameter)  passing 
through  the  measurement  volume  formed  by  the  Intersection  of  the  laser  beams 
in  Figure  3.1  and  to  determine  the  particle  size  by  statistical  analysis  of 
the  signal.  Although  the  system  shown  is  physically  similar  to  the  familiar 
laser  Doppler  veloclmeter,  important  differences  in  its  implementation  form 
the  basis  of  the  present  approach.  As  explained  in  Chapter  2,  the  Doppler 
signal  obtained  from  the  scattered  light  only  contains  time  (but  not  velocity) 
information  which  must  be  analyzed  in  a  statistical  sense.  This  concept  is 
that  the  mean  excursion  distance  for  the  particles  in  Brownian  motion  is  much 
smaller  than  the  wavelength  of  visible  light,  as  shown  in  Figure  2.4.  As 
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discussed  In  greater  detail  In  Chapter  2,  the  consequence  of  this 
characteristic  of  the  measurement  Is  that  the  magnitude  of  signal  variations 
is  not  meaningful  (because  the  motion  may  be  occurring  in  different  regions  of 
the  LDV  fringe  pattern).  The  frequency  characteristics  of  the  signal  are 
significant,  however,  in  that  they  can  be  related  to  the  statistics  of 
particle  motion.  The  identification  of  such  relationships  is  the  motivation 
for  the  computer  simulation  of  the  signal. 

The  basis  of  the  computer  simulation  is  the  probabilistic  solution  (Eq. 
2.8)  of  the  Langevln  equation.  Using  the  probability  distribution  of  Eq. 
(2.8),  we  may  generate,  using  a  pseudo-random  number  generator,  the  velocity 
history  of  w  with  the  correct  statistical  properties.  In  general,  if  R  is  a 
random  number  with  uniform  distribution  in  the  interval  [0,1],  and  p(x)  is  the 
probability  distribution  of  a  variable  x,  R  may  be  mapped  into  x  by  equating 
the  probability  of  finding  x  in  [xj ,  x2], 
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In  particular  let  xj  -  -®  and  R*  ■  0;  then 
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For  the  distribution  of  Eq.  (2.8),  it  is  convenient  to  rewrite  the  expression 
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Then  Eq.  (10)  becomes 
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from  which  wn  may  be  solved  from  the  random  number  R.  The  solution  of 
Eq.  (3.6)  is 


where 


£  -  erf"1  (2R-1) 

n  n 


(3.7) 
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To  summarize,  the  algorithm  for  generating  the  Brownian  velocity  history 
works  as  follows.  The  Initial  velocity  wj  is  obtained  randomly  from  the 
Maxwellian  distribution 
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The  subsequent  velocities  wQ,  with  sampling  time  At,  are  then  given  by 
Eqs.  (3.7)  and  (3.8). 


3.2  The  Doppler  Signal 


Consider  a  crossed  laser  beam  system  with  beam  waist  rQ  as  shown  in 
Figure  3.1.  The  intensities  of  the  two  beams  scattered  by  the  particle  at  x  - 
(x,y,z)  (in  a  coordinate  system  centered  on  the  probe  volume)  are 
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Here  the  unit  propagation  vectors  are 
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Therefore  Che  expressions  for  1^  and  I2  become 
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(3.14,  3.15) 

The  mean  particle  trajectory  Is 


x(t)  -  x  +  v  t  ,  (3.16) 

—  — o  — 

where  Xq  is  the  reference  point  and  v  Is  the  mean  velocity.  Substitution  of 
x(t)  ■  (x(t),  y(t),  z(t))  from  Eq.  (3.16)  Into  Eqs.  (3.11)  and  (3.12)  would 
give  the  scattered  Intensity  profile  as  the  particle  traverses  the  beam 
volume. 

At  the  detector,  because  of  the  interference  between  the  two  beams,  the 
light  collected  Is  a  modulated  Doppler  signal.  The  signal  may  be  written  as 

t 

S  -  I.  +  I,  +  2  /I.I_  cos  Uk  •  /  vdt+*)  ,  (3.17) 
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where 

0 

Akz  -  yL  (2  sin  •—)  ,  (3.18) 

and  the  z  velocity  component  contains  both  a  mean  and  random  velocity, 

v  -  v  +  w(t)  .  (3.19) 

z  z 

The  phase  shift  <p  is  a  constant  depending  on  the  reference  point  definition. 
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The  power  received  at  the  detector  nay  therefore  be  written  as 
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(fj  +  f2  +  2/f jf 2  cos  •  2  sin  ^  •  [v(  t  t  /  w(t)dt])} 


(3.20) 


The  first  term  represents  the  absolute  Intensity  function  which  is 
proportional  to  the  beam  Intensity  I0,  the  aperture  and  transmissivity 
function  H  of  the  receiving  optics,  and  the  scattering  cross-section 
Integrated  over  the  appropriate  field  of  view.  The  functions  fj  and  f2  are 
defined  by  the  right  hand  sides  of  Eqs.  (3.11)  and  (3.12).  The  cosine  term 
represents  the  heterodyned  Doppler  signal. 


3.3  Simulation  Results  and  Discussion 

The  Brownian  velocity  of  a  0.1  urn  diameter  spherical  particle  with  a 
specific  gravity  of  2  In  300  K  air  at  atmospheric  pressure  has  been  simulated 
using  the  Monte-Carlo  calculation  described  above.  The  velocity  history  Is 
shown  In  Figure  3.2  for  500  velocity  samples  with  a  sampling  time  of  10  ns 
each.  The  velocity  statistics  are  shown  In  the  histogram  In  Figure  3.3.  The 
result  agrees  with  the  Maxwellian  distribution  which  Is  also  plotted  In  the 
figure. 

Calculations  were  then  carried  out  to  examine  the  behavior  of  the  Doppler 
signal,  which  Is  described  by  Eq.  (3.20).  The  term  of  interest  is 


s  -  |{fl  +  f2  +  coe  (f1  •  2  sin  ^  [vzt  +  /  w(t)dt]  +  <(.)}  ,  (3.21) 
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Figure  3.2.  Simulated  Brownian  Velocity  History  of  a  0 
Particle  (Specific  Gravity  ■  2)  in  Air  at  ' 
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with  f1  and  f2  defined  by  the  right  hand  sides  of  Eqs.  (3.11)  and  (3.12).  To 
obtain  the  most  favorable  signal,  the  particle  is  assumed  to  have  mean 
velocity  in  the  y  direction  only,  and  is  traversing  the  center  of  the  sample 
volume  along  the  y  axis.  The  reference  phase  value  Is  set  arbitrarily  to 
(-0.5* )  which  corresponds  to  eligning  the  y  axis  along  the  "boundary”  of  the 
bright  and  dark  fringes  in  the  bean  sample  volume.  The  laser  beam  waist  r0 
is  50  um  and  the  beam  crossing  angle  is  set  at  64  -  10°.  These  values  are 
typical  of  a  practical  system. 

Under  the  above  conditions,  the  simulated  Doppler  signal  for  a  0.1  um 
particle  with  vy  “  1  m»s_1  traversing  the  center  part  of  the  beam  volume, 
starting  from  x  ■  (0.,  -10  um(  0.),  is  shown  in  Figure  3.4(a).  The  horizontal 
axis  Is  the  number  of  time  steps,  with  each  time  step  equal  to  10  ns.  The 
signal  fluctuation  due  to  the  Brownian  motion  is  clearly  detectable. 

The  simulation  is  repeated,  under  the  same  conditions,  with  a  0.01  um 
diameter  particle.  The  simulated  signal  for  the  smaller  particle  Is  shown  In 
Figure  3.4(b).  The  fluctuation  Is  of  a  much  higher  level  than  those 
exhibited  by  the  larger  diameter  particle.  Furthermore,  because  of  the 
behavior  of  the  cosine  function,  the  modulation  Is  reduced  when  the  argument 
of  the  cosine  Is  close  to  +  n*,  n  -  0,1,2,....  For  the  signal  In  Figure 
3.4(b),  the  pertlcle  Is  very  close  to  the  center  of  the  beam  volume  as  the 
absolute  Intensity  Is  approximately  constant.  Therefore,  the  signal  has  a 
magnitude  of 
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1  +  cos  •  2  sin  ~  [v  t  +  /  w(t)dt]  + 


(3.22) 


Thus,  at  regions  where  S  *1  or  2,  the  magnitude  of  the  modulation  Is 
reduced.  This  effect  is  evident  In  Figure  3.4(b). 


The  root-mean-square  value  <z>nBS  of  the  Brownian  displacement  Is  shown  In 
Figure  3.5.  For  particles  of  100  to  1000  A  diameter,  <zm>/X  ~  0.1. 

Therefore  the  nns  displacement  is  much  less  than  the  fringe  spacing  k'“  .  As 
a  result,  the  length  scale  A  does  not  play  any  role  in  the  determination  of 
the  Brownian  velocity  w.  The  signal,  therefore,  only  contains  Information 
about  the  time  scale  and  not  any  length  scale,  and  the  determination  of  a 
meaningful  velocity  Is,  in  principle,  not  possible. 

The  essence  of  the  above  discussion  is  illustrated  in  Figure  3.6,  which 
exhibits  the  transfer  function  of  the  detection  system  [the  right  hand  size  of 
Gq.  (3.31)].  For  typical  Brownian  motion  modulation  (AB),  the  signal 
modulation  is  (A'B').  The  signal  modulation  is  drastically  reduced  where  the 
slope  of  the  transfer  function  vanishes  (CD).  The  amplitude  of  the  modulation 
Is  a  function  of  the  position  of  the  particle  relative  to  the  fringe  pattern. 
Since  this  is  not  known  a  priori,  the  amplitude  modulation,  therefore,  does 
not  provide  useful  information  on  the  nature  of  the  Brownian  motion. 
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Figure  3.6.  Modulation  of  the  Doppler  Signal  as  a  Function  of 
Particle  Position 


When  there  is  a  fluid  velocity  component  vz  perpendicular  to  the  fringe 
pattern,  the  particle  position  will  be  represented  by  the  trajectory  (EF), 
with  the  Brownian  modulation  such  as  (AB)  superimposed  on  top.  The  signal 
will  exhibit  full  modulation  (E'F')  (regular  Doppler  signal)  with  the  Brownian 
motion  modulation  (A'B')  superimposed  on  top.  When  v2  is  large  compared  to 
the  Brownian  velocity  w,  it  may  be  difficult  to  recover  the  Brownian  motion 
modulation  from  the  regular  Doppler  signal.  This  effect  is  illustrated  in 
Figures  3.7(a)  and  3.7(b)  in  which  velocity  vz  of  0.1  and  0.3  m«s-1 
respectively  are  superimposed  onto  the  Brownian  velocity.  In  the  latter  case, 
the  Brownian  motion  signal  is  almost  completely  masked  and  the  signal 
resembles  the  usual  Doppler  velocimeter  signal.  For  small  particles  with  a 
higher  Brownian  velocity,  good  quality  signal  may  still  be  obtained  with  a 
significant  vz.  This  is  illustrated  in  Figure  3.8  by  the  simulation  of  the 
signal  from  a  0.01  p.m  particle  with  vz  *  0.3  m^s-1  .  The  mean  Brownian 
velocity  w,  where 


(3.27) 


is  approximately  0.1  m* s-  for  a  0.1  ^m  diameter  particle  at  a  gas  temperature 
of  300  K.  Thus  one  expects  that  for  particles  in  the  0.01  to  0.1  ^  diameter 
range,  reasonable  signal  quality  may  be  obtained  if  "vz  is  less  than 
0.1  m*  s-1  . 

The  above  discussion  indicates  that  the  strategy  for  analyzing 
characteristics  of  the  Brownian  motion  lies  in  analysis  of  the  time 
information.  The  fundamental  time  scale  of  the  fluctuation  is  the  relaxation 
time  1/3*  The  information  on  the  particle  size  .uay  be  calculated  from  3  using 
Eqs.  (2.4)  and  (2.5).  Furthermore,  the  signal  must  be  filtered  to  remove  the 
regular  Doppler  component  (vz)  for  proper  interpretation. 


SIGNAL 


THE  step 

(10  ns/steo) 

(a) 


Tir-E  STEP 
(10  ns/step) 


(5) 

Figure  3.7.  Simulated  Brownian  Motion  Signals  for  a  0.1  um  Diameter  Particle 
with  Fluid  Velocity  Component  vz  Superimposed 

(a)  vz  =•  0.1  m  s-1 

(b)  vz  -  0.3  m  s-1 


3 .4  Processing  of  the  Brownian  Motion  Sensor  Signal 

As  discussed  above,  only  time  information  may  be  obtained  from  the 
Brownian  motion  sensor.  Since  the  time  scale  relevant  to  the  Brownian  motion 
is  much  faster  than  any  other  time  scales  (e.g.,  the  transit  time  of  the 
particle  due  to  its  mean  velocity),  the  latter  may  be  filtered  out  easily.  In 
the  following  discussion  we  shall  assume  that  the  signal  has  already  been 
filtered  to  contain  the  Brownian  information  only.  The  appropriate  statistics 
of  the  temporal  information  signal  are  discussed. 

As  depicted  in  Figure  3.6,  the  signal  amplitude  is  a  function  of  the 
absolute  position  of  the  particle  within  the  probe  beam  volume.  Therefore, 
the  relevant  information  is  contained  in  the  time  between  the  extrema  of  the 
signal.  At  an  extremum,  there  is  a  change  of  direction  of  the  particle,  and, 
therefore,  the  extremum  represents  a  zero  crossing  of  the  particle  velocity. 
Thus,  the  Brownian  motion  may  be  analyzed  in  terms  of  the  time  between  extrema 
of  the  signal  (zero  crossing  of  the  first  derivative).  The  particular  method 
of  implementation  depends  on  the  hardware  configuration.  Nevertheless  the 
principle  remains  the  same.  Also,  the  analysis  considers  only  the  mean  time 
between  extrema,  and  does  not  to  attempt  analyze  the  full  spectral  content  of 
the  signal,  so  that  possible  aliasing  is  not  an  issue.  (The  method  amounts  to 
doing  statistical  sampling  rather  than  fixed  interval  sampling,  so  on  the 
average,  the  higher  frequency  content  is  "represented.") 

The  signal  processing  is  simulated  on  the  computer  using  ‘he  Monte  Carlo 
technique  described  previously.  The  Brownian  motion  of  a  particle  with 
specific  gravity  of  two,  embedded  in  room  temperature  air,  is  analyzed.  The 
mean  time  between  extrema  "t  as  a  function  of  the  particle  diameter  d  is  shown 
In  Figure  3.9.  The  values  of  t  do  not  seem  to  follow  a  simple  linear 
relationship  with  d,  yet  the  only  time  scale  involved  is  the  relaxation  time 
x  relax  *  1/f3*  With  0  proportional  to  1/d  in  the  free  flow  regime,  the 
relaxation  time  should  be  proportional  to  d.  Furthermore,  depending  on  the 
sampling  interval  At,  the  mean  time  t  is  different.  The  relevant  parameter 
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should  therefore  clearly  involve  At.  This  issue  is  resolved  when  the  data  are 
replotted  against  /At/(3  (with  the  dimension  of  time)  (Figure  3.10).  The 
simulation  points  collapse  on  to  a  single  straight  line  passing  through  the 
origin.  The  relationship  between  sampling  and  adequate  statistics  is  also 
illustrated,  at  large  values  of  /At/g ,  a  long  string  of  data  is  required  to 
produce  a  good  estimate  of  t.  For  small  values  of  /At/B ,  the  sampling 
interval  of  At  *  10  ns  is  not  adequate  to  resolve  the  Brownian  motion  time 
scale,  and  as  a  result,  t  is  overestimated. 
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Figure  3-10.  Correlation  Between  the  Mean  Time  Between  Signal  Extrema 
and  the  Particle  Relaxation  Time,  6_1 


4.  THEORETICAL  ANALYSIS  OF  THE  DOPPLER  SIGNAL 


In  Che  last  section,  a  Monte  Carlo  simulation  of  the  signal  from  the 
Brownian  motion  sensor  was  carried  out.  The  result  of  interest,  namely,  the 
time  between  zero  crossings  of  the  Brownian  velocity,  was  found  to  be 
dependent  on  the  step  size  At  of  the  simulation.  This  dependence  is  actually 
not  surprising  because  At  represents  essentially  the  bandwidth  of  the  actual 
signal  detection  system,  which,  being  a  practical  device,  should  have  a  finite 
bandwidth.  Since  the  Brownian  motion  is  a  statistical  process,  the  behavior 
of  any  signal  derived  from  the  Brownian  motion  is  statistical,  and  the 
spectral  properties  of  the  signal  will  depend  on  the  bandwidth  of  the  signal 
detection  system.  These  ideas  are  further  explored  in  this  chapter.  Also  an 
analytical  approach  is  attempted  in  order  to  provide  a  theoretical  basis  for 
the  signal  processing  method  used  in  Chapter  3. 


The  Statistical  Nature  of  the  Brownian  Velocity 


The  Brownian  velocity  of  a  particle  in  thermal  equilibrium  in  a 
background  fluid  is  described  by  the  Langevin  equation  (see  Chapter  2).  For 
the  observed  velocity  component  u  (in  the  direction  of  Ak),  the  Langevin 
equation  is 
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where  is  a  relaxation  time  governed  by  the  mass  m  and  the  mobility  e  of 
the  particle  in  the  fluid. 
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The  statistical  nature  of  u  Is  due  to  the  random  driving  force  f,  which  is  the 
result  of  the  impulse  arising  from  each  collision  of  the  particle  and  a 
molecule  of  the  background  fluid.  Because  these  collisions  occur  randomly,  f 
may  be  considered  to  be  a  white  noise  function  with  a  constant  spectral 
density  Sf, 


<  f ( t ) f ( t+r )  >  -  Sf6(t) 


(4.3) 
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The  time  scale  associated  with  the  fluctuation  of  f  is  of  the  order  of 
the  time  between  the  molecular  collisions  rc,  which,  for  a  0.1  pun  particle 
in  a  gas  medium  at  STP  is  approximately  10“14  s.  The  time  response  of 
practical  Instruments  used  to  observe  the  Brownian  motion  velocity  is  usually 
much  slower  than  this  collision  time.  The  observations,  are,  therefore,  the 
result  of  averaging  u  over  a  response  time  scale  At  associated  with  the 
instrument,  with  At  >  xc.  The  observed  value  of  u  is  then  a  random 
variable,  and,  since  it  is  the  aggregate  of  many  similar  events,  it  has  a 
Gaussian  distribution,  by  the  central  limit  theorem.  If  the  particle  has  a 
velocity  uQ  at  one  instant  of  time,  the  conditional  probability  distribution 
of  u  at  a  time  At  later,  (so  that  the  value  of  u  is  the  aggregate  sum  of  all 
the  collisions  in  At),  is  given  by [4.1] 
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Note  that  the  distribution  is  parameterized  by  the  "memory"  factor  h 


(4.6) 


4-2 


H 


r  O  C  >■  s'  •.*  V  s*  s'  s'  s'  s'  s'  s*  s”  s'  s'  s'  s'  s'  v*  s'  s'  s'  s'  s'  s'  s'  s'  s'  ■’.‘'.'T.-'.  ,  -■ 


*  f  f-crjw  U--SS  .V, 


and  the  width  a.  The  observed  velocity  u(t)  may  be  viewed  as  observations  at 
discrete  points  in  time,  tn  “  nAt.  When  the  time  step  At  is  large  compared 
to  the  relaxation  time  (pAt  >1),  the  velocities  U£  are  independent  of  each 
other,  and  each  obeys  a  Maxwellian  distribution  at  temperature  T.  In  the 
Brownian  motion  sensor  described  here,  since  the  velocity  amplitude  is  not 
measured,  all  the  information  is  contained  in  the  measurement  of  the  value  of 
p.  Therefore,  the  instrument  response  time  At  should  be  such  that  PAt  <  I. 
After  p  is  determined,  the  particle  physical  properties  (e,m)  may  be 
determined  from  Eq.  (4.2). 


4.2  The  Distribution  of  the  Zeros  of  the  Velocity 


The  Brownian  motion  sensor  measures  the  time  between  successive  zeros  of 


the  Brownian  velocity  u  (in  the  direction  of  Ak) .  The  time  history  of  the 
velocity  values,  u(t),  as  observed  by  an  instrument  with  an  integration  time 
At,  may  be  thought  of  as  a  time  series  un  at  time  nAt.  From  Eq.  (4.4)  the 
values  of  each  un  are  governed  by  the  probability  distribution 
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where  the  un_j  appears  as  a  parameter.  We  want  to  determine  the  probability 
P(N)  of  obtaining  a  time  between  two  successive  zeros  of  NAt. 


The  problem  of  determining  the  distribution  of  zeros  of  a  random  variable 
has  been  studied  extensively  in  the  context  of  shot  noise  in  electronic 
devices,  and  in  queueing  theory.  The  classic  review  of  the  subject  is  given 
in  Reference  4.2.  We  shall  derive  the  expression  for  P(N)  Independently 
here.  To  connect  our  result  to  that  in  the  literature  of  noise  current,  the 
random  driving  force  f  in  Eq.  (5)  may  be  considered  as  the  randomly  arriving 
electron,  and  the  observed  "current”  u  is  the  result  of  two  transfer  functions 
(or  filters)  in  series.  The  first  transfer  function  is  that  associated  with 


the  Langevin  equation.  Formally 


VO  '.A  _  .W 


V.H  V'  S';  V 


u(t) 


/ 

—00 


(p+ioj)  2  it 


(4.8) 


where  f  is  the  Fourier  Transform  of  f.  Physically,  this  transfer  function  is 
due  to  the  averaging  of  the  effect  of  f  by  the  damping  of  the  particle  in  the 
fluid.  The  second  transfer  function  is  due  to  the  averaging  by  the 
instrument.  Therefore  the  overall  transfer  function  is 
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It  should  be  noted  that  no  closed  form  solutions  for  P(N)  have  been  worked 
out,  except  for  rather  trivial  cases.  The  difficulty  is  attributed  [4.2]  to 
the  lack  of  an  analytical  solution  to  the  N  dimensional  integral  involved  (see 
Eq.  (4.10). 

We  shall  derive  the  expression  for  P(N)  here.  The  process  of  obtaining  a 
successive  zero  in  N  steps  may  be  described  as  the  following.  At  time  zero 
(n  *  0),  the  particle  has  a  very  small  velocity  6u,  which,  within  the 
amplitude  resolution  of  the  instrument,  may  be  considered  as  zero.  In  the 
following  argument,  6u  is  taken  as  positive.  The  argument  will  be  exactly  the 
same  if  6u  is  negative.  We  shall  define  a  zero-crossing  event  at  t  •  NAt, 

(n  *  N),  as  an  event  with  velocity  greater  than  zero  at  t  *  (N  -  1)  At, 

(n  ■  N  -  1),  and  velocity  less  than  or  equal  to  zero  at  t  »  NAt.  The  time 
resolution  of  this  definition  of  zero-crossing  is  At,  and  so  the  meaning  of 
such  an  event  for  N  ■  1  is  not  precise.  Then,  the  condition  for  occurrence  of 
a  successive  zero  in  N  time  steps  is  to  have  the  velocities  be  positive  for 


n  -  1,  ,  N-l ;  and  the  velocity  be  negative  or  equal  to  zero  for  n*N. 

Therefore 


P(N)  -  /  du^  p  (Uj,  0;  At)  /  du2  p  (u2>  u^  At)  •••  (4.10) 
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Using  the  formula  for  p  in  Eq.  (4.7),  and  changing  the  variable  from  u  to 
£  -  u/o,  the  above  probability  becomes 
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(4.11) 


Note  that  a  does  not  appear  in  this  expression.  This  is  because  the 
zero-crossing  time  pertains  to  the  zeros  of  the  velocity,  and  is  therefore 
independent  of  the  velocity  scale.  (The  value  a  is  a  velocity  scale  measuring 
the  width  of  the  velocity  distribution). 

If  the  integration  time  At  is  large  (0At  »  1),  h  will  approach  zero  and 
Eq.  (4.11)  reduces  to  a  binomial  process 

1  N 

P(N)  -  (j)  (4.12) 


The  average  value  of  N  is 
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In  this  case,  the  data  values  corresponding  to  the  successive  zeros  of  the 
velocity  become  a  completely  random  sequence,  and  do  not  contain  information 


about  h. 


For  small  values  of  h,  the  multiple  integral  may  be  worked  out  as 
follows:  An  operator  In  is  defined  as 


*n  <f(V> 
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The  quantity  F,n_i  enters  as  a  parameter  in  the  above  definition.  Then 


\  +  +  0(h2> 


(4.15) 


(5„h)  *  I  J  +  0<h2) 


(4.16) 


The  expression  for  P(N),  Eq.  (4.11),  after  evaluating  the  last  integral  on  the 
right  hand  side,  becomes 


P(N)  41!  I2  •••  Vi  u  +  J  sN_ih)  +  0(h2) 


(4.17) 


The  expression  (4.15)  and  (4.16)  may  be  used  repeatedly  in  (4.17).  The  first 
few  expressions  are: 
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The  above  expression  is  valid  for  N  ^  2.  For  N  »  1 ,  since  the  velocity  at 
n  =  0  is  preconditioned  to  be  «  0, 
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Compared  to  the  binomial  distribution  (P(N)  *  2“N),  which  is  the  limit 
for  h  ■  0,  P(N)  decreases  with  h  for  N  <  3  and  increases  with  h  for  N  >  3. 

This  result  is  illustrated  in  Figure  4.1.  Since  the  velocity  at  n  •  0  is 
zero,  the  velocity  u^  at  n  ■  1  is  sampled  from  a  Gaussian  centered  at  zero. 

The  value  of  P(N  »  1)  is  the  probability  that  this  velocity  is  negative,  and 
P(l)  is  therefore  1/2.  Say  the  velocity  at  n  •  1  is  Uj(>  0),  then  the 
velocity  u2  at  n  ■  2  is  sampled  from  a  Gaussian  centered  at  u^h.  The 
probability  of  u2  being  negative  is  the  shaded  area  in  Figure  4.1  under  the 
Gaussian  at  n  *  2,  and  this  value  is  <  1/2  for  h  >  0.  Therefore  P(2) 

<  (1/2)2.  For  N  *  3,  the  probability  of  obtaining  a  consecutive  ;ero  here  is 
represented  by  the  shaded  area  under  the  Gaussian  centered  at  u2h  at  n  ■  3 
(Figure  4.1).  This  area  has  to  be  summed  over  all  the  realizations  of  the 
possible  values  of  positive  u2 ,  which  is  sampled  from  the  Gaussian  centered  at 
Ujh,  which  again  has  to  be  summed  over  all  the  realizations  of  u  ^ .  Although 
the  shaded  area  at  n  ■  3  is  less  than  1/2,  there  are  "more"  realizations  than 
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the  h  =  0  case  because  u2  is  sampled  from  the  shifted  Gaussian.  These  two 
effects  balance  each  other  to  0(h2),  and  P(3)  is  equal  to  (1/2)3  +  0(h2). 

For  N  >  3,  P(N)  is  given  by  Eq.  (22),  and  is  greater  than  (1/2)N.  The 
values  of  P(N)  are  plotted  versus  N  in  Figures  4.2  and  4.3. 

It  is  not  possible  to  obtain  a  closed  form  solution  of  P(N)  for  arbitrary 
h  (h  <  0  <  1  corresponding  to  0  <  6At  <  <*>)  because  of  the  difficulty  in 
evaluating  the  multiple  integral  in  Eq.  (4.7).  Direct  numerical  calculation 
the  Integral  is  not  feasible  for  any  large  N,  say  N  >  5.  Monte  Carlo 
evaluation  of  the  integral  would  be  equivalent  to  our  earlier  simulation  of 
the  Brownian  motion  in  Chapter  3.  The  results  of  Chapter  3  will  be  reexamined 
in  this  context  here. 

The  general  behavior  of  the  mean  number  of  steps  N  between  consecutive 
zero-crossings  of  the  Brownian  velocity  is  shown  in  Figure  4.4.  For  h  *  0,  N 
is  2  (from  Eq.  (4.13)).  For  small  h, 


l  N  (i)  (1  +.  (N-3)  — )  +  0(h2) 
N=1  Z 


l  <1  -?)  +  !  (tj)  (?)  +  °(h2> 

N-l  2  N-l  2 


N  -  0(h2)  (4.20) 

In  the  last  expression,  we  made  use  of  J(N/2N)  «  2  and  £(N2/2N)  *  6. 

Therefore  N  *  0(h2)  for  small  h.  For  large  relaxation  time  (6  +  0,  8At  +  0, 
and  h  +  1),  the  particle,  once  it  acquires  a  positive  (or  negative)  velocity, 
will  retain  that  velocity  for  a  large  time.  Consequently  N  +  =°. 


The  results  of  the  Monte  Carlo  simulation  in  Chapter  3  have  been 
replotted  in  the  context  of  the  present  analysis  in  Figure  4.5,  and  confirm 
the  theory.  The  results  in  Figures  4.4  and  4.5  show  that  to  obtain  h,  (and 
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Figure  4.4.  General  Behavior  of  the  Expectation  Value  of  the  Number 
of  Time  Steps  Between  Consecutive  Zero-Crossings  of  the 
Brownian  Velocity,  as  a  Function  of  the  Memory  Parameter 
h  ■  exp  <-84t). 
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Figure  4.5.  Monte  Carlo  Simulation  Result  of  the  Mean  Time  Step  (N) 

Between  Zero-Crossings  of  the  Velocity,  (a)  as  a  Function 
of  64t;  (b)  as  a  Function  of  h  - 


therefore  the  values  of  B,  which  would  lead  to  the  values  for  £/m) ,  it  Is 
necessary  to  have  a  time  resolution  of  At  such  that  BAt  <  1  (or  h  >  e" ;'at 
close  to  1) . 


4.3  Summary 

A  theoretical  analysis  of  the  probability  of  obtaining  a  time  t  between 
consecutive  zero-crossings  of  the  Brownian  velocity  u  has  been  presented. 

Since  u  is  a  random  variable,  the  probability  P(t)  is  necessarily  a  function 
of  the  integration  time  At  of  the  Instrument  used  to  observe  u.  The  analysis 
has  been  carried  out  in  terms  of  discretizing  the  time  in  steps  of  At, 
tn  ■  nAt,  and  the  resulting  probability  P(N),  (N  -  T/ At),  has  been 
derived.  The  values  of  P(N)  are  found  to  be  dependent  on  the  "memory"  h  (h  = 
exp  (-BAt))  of  the  Brownian  velocity,  where  6-1  is  the  relaxation  time.  For 
very  short  memory  (h  <  1),  the  probability  P(N)  tends  to  the  binomial  limit 
P(N)  -  1/2n.  The  expectation  value  N  (  =  J  NP(N)  assumes  the  value  of  2  in 
the  binomial  limit  and  2  +  0(h2)  for  small  h.  For  large  relaxation  time  (B  ♦ 
0) ,  N  *  <*>.  To  obtain  accurate  information  from  the  measurement  of  N  about  h, 
and  therefore,  the  value  of  S  which  is  related  to  the  physical  properties  of 
the  particle  suspended  in  the  fluid,  it  is  necessary  to  have  BAt  <  1.  Since 
the  result  is  a  function  of  the  Integration  time  At,  which  is  not  a  precisely 
defined  quantity  for  the  measuring  Instrument,  the  Brownian  motion  sensor  must 
be  a  calibrated  Instrument. 


References  for  Chapter  4 


1.  Chandrascklar,  S.,  "Stochastic  Problems  in  Physics  and  Astronomy, 
Reviews  of  Modern  Physics,  Vol.  15,  1,  1943. 

2.  Rice,  S.O.,  "Mathematical  Analysis  of  Random  Noise,”  Bell  System 
Technical  Journal,  Vol.  23  and  24. 


5.  LABORATORY  EXPERIMENT  FOR  BROWNIAN  MOTION  DETECTION 


In  this  chapter,  the  signal  to  noise  requirement  for  detecting  the 
Brownian  motion  of  a  small  particle  using  an  interferometric  system  is  first 
estimated.  Then  the  experimental  set  up  and  the  data  acquisition  systems  are 
described.  Finally,  the  signal  sampled  from  the  detection  system  is 
presented. 

5.1  Signal  to  Noise  Requirement  of  Brownian  Motion  Detection 

The  signal  from  the  photomultiplier  has  a  noise  component  due  to  the 
random  arrival  of  the  photo-electrons  (shot  noise),  a  deterministic  frequency 
component  due  to  the  convection  of  the  particle  through  the  fringed  volume 
formed  by  the  intersecting  laser  beams  (regular  Doppler  signal),  and  a 
statistical  component  due  to  the  Brownian  motion  of  the  particle.  Since  the 
displacement  of  the  particle  due  to  Brownian  motion  is  small  compared  to  the 
fringe  spacing  (see  Chapter  2),  all  the  information  that  can  be  obtained  from 
the  signal  on  the  Brownian  motion  characteristics  is  contained  in  the  small 
fluctuation  which  is  superimposed  on  top  of  the  regular  Doppler  signal.  To  be 
able  to  detect  this  information,  the  signal  fluctuation  due  to  Brownian  motion 
must  be  much  larger  than  that  due  to  shot  noise  within  the  bandwidth  of  the 
detector. 

The  bandwidth  Af  of  the  detector  should  be  able  to  resolve  the  relaxation 
time  xreiax  of  the  particle.  For  particles  0.02  to  0.2  ^m  in  diameter, 

Trelax  x  to  3  x  10-7  s.  Therefore  Af  is  >  3  x  107  Hz.  The  shot 


noise  current  is 


For  a  heterodyne  system,  the  statistics  limited  signal  to  noise  ratio  is 
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(5.2) 


where  n  is  the  number  of  photo-electrons  collected  per  unit  time  and  Af  is  the 
bandwidth.  The  sine  factor  is  due  to  the  fact  that  the  excursion  of  the 
particle  due  to  Brownian  motion  is  small  compared  to  the  fringe  spacing  \f, 
hence  the  modulation  of  the  heterodyned  signal  is  not  complete,  and  the 
modulation  depth  is  proportional  to  the  maximum  excursion  distance  X,,, 
occurring  during  xreiax«  Furthermore,  the  bandwidth  should  be  fast  enough 
to  capture  Tre^ax  for  the  smallest  particle  of  interest.  The  relaxation 
time  for  a  0.02  ^m  particle  is  4  x  10“8  s.  Therefore  a  bandwidth  of  25  MHz  is 
sufficient.  The  data  acquisition  system  to  be  described  has  a  bandwidth  of  30 
MHz,  and  this  figure  is  used  for  the  signal  to  noise  calculation.  (Note  that 
to  resolve  treiax,  the  sampling  rate  should  be  at  least  twice  the  bandwidth 
according  to  the  Nyquist  theorem.  In  our  data  acquisition  system,  the 
sampling  rate  was  100  MHz). 

In  terms  of  the  parameters  of  the  optical  system,  the  signal  to  noise 
ratio  for  a  collection  optical  system  oriented  at  a  90°  scattering  angle  may 
be  expressed  as 


S/N 
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sin(  2xxm/\^) 
Af 


(5.3) 


where  I  is  the  incident  beam  power  (photon/sec),  is  beam  cross  section 
area,  xr  is  the  collection  optics  efficiency,  r\  is  the  quantum  efficiency, 
dQ  is  the  collection  solid  angle,  and  dc/dQ  is  the  differential  scattering 
cross  section.  For  Rayleigh  scattering 
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The  dJ  size  dependence  is  due  to  the  variation  of  the  Rayleigh  scattering 
cross  section  (~d6). 

For  particles  in  the  0.02  to  0.2  (±m  range,  Eq.  (5.5),  which  is  based  on 
Rayleigh  scattering,  is  not  strictly  valid.  Comparisons  with  the  exact  Mie 
calculations  at  0.2  jam  for  a  laser  wavelength  of  0.5  um,  however,  showed  that 
the  scattering  amplitudes  calculated  by  the  Rayleigh  formula  are  less  than  30% 
off.  For  an  order  of  magnitude  S/N  estimate  in  the  0.02  to  0.2  p.m  range, 
therefore,  Eq.  (5.5)  is  sufficient.  The  signal  to  noise  ratios  for  systems 
with  laser  power  density  of  1010  W/m2  and  109  W/ra2  for  each  of  the 
Intersecting  beams,  and  with  f / 1.3  collecting  optics,  xr  *  0.3  and  q  *  0.2 
are  shown  in  Figure  5.1.  The  refractive  indexes  used  in  the  calculations 
were  1 .5  +  Oi,  which  corresponds  to  that  of  latex  spheres,  and  4.5  -  1.9  1, 
which  corresponds  to  solid  carbon  particles.  (The  latter  value  is  used  as  a 
simulation  for  the  S/N  in  soot  measurements.) 

For  measurements  on  latex  particle  used  as  a  calibration  procedure,  the 
estimated  signal  to  noise  is  marginal,  ~1  to  10,  even  for  particles  up  to 
0.2  (im.  The  actual  S/N  attainable,  therefore,  is  very  much  dependent  on  the 
actual  apparatus. 


>aratus  Development 


The  observation  of  the  Brownian  motion  of  individual  small  particles 
poses  several  experimental  design  challenges  due  to  the  small  scattering 
cross  section  of  the  submicron  particles  combined  with  the  short  observation 
time  required  to  resolve  the  Brownian  motion.  As  a  result,  the  experiment  was 
severely  limited  by  signal  to  noise  considerations,  and  much  effort  focused  on 
refining  the  optical  system  to  improve  the  signal  to  noise  ratio.  Several 
versions  of  the  apparatus  were  developed  in  an  evolutionary  manner. 

The  optical  system  which  had  been  first  assembled  for  this  study  is  shown 
schematically  in  Figure  5.2.  A  Spectra-Physics  Model  165  argon  ion  laser  with 
a  nominal  single-line  power  of  approximately  2.5  W  which  yields  a  true 
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Figure  5.1.  Estimated  Signal  to  Noise  Ratio  of  Brownian  Notion  Detection 
System.  (The  calculation  is  not  valid  for  d_>  0.2  jm  because 
the  Rayleigh  scattering  formula  was  used  for  the  scattering 
cross-section.) 


TEMqo  output  with  optimum  focusing  characteristics  was  used.  To  obtain  a 
clear  fringe  pattern  at  the  measurement  volume,  it  was  necessary  to  insert 
a  Spectra-Physics  Model  589  air  spaced  etalon  into  the  laser  cavity  to 
restrict  laser  oscillation  to  a  single  longitudinal  mode.  This  provided  the 
stable  mode  amplitude  and  increased  coherence  length  required  in  our  optical 
measurement  system,  but  at  the  expense  of  a  75%  decrease  in  laser  power.  The 
vertically  polarized  output  of  the  laser  was  expanded  by  25x  and  collimated  to 
a  beam  size  of  approximately  3  cm.  The  beam  expansion  lead  to  tighter 
focusing  and  a  resultant  higher  power  density  at  the  sample  volume  than  could 
be  achieved  with  a  smaller  diameter  beam.  An  adjustable  iris  diaphragm  was 
used  after  the  beam  expander/collimator  to  select  the  beam  diameter. 

The  expanded  beam  was  split  into  two  beams  (b^  and  hjO  by  a  50/50 
cube  beamsplitter,  which  reflected  50%  of  the  incident  light  and  transmited 
35%.  An  increasing  intensity  difference  in  two  beams  when  focused  to  a  point 
to  form  interference  fringes  resulted  in  a  decreasing  clarity  of  the  fringe 
pattern,  which  was  detrimental  to  the  signal  to  noise  ratio  in  our  optical 
system.  Therefore,  three  reflectors  were  used  to  decrease  the  intensity  of 
the  reflected  beam  (=»  10%  light  loss  off  of  each  mirror)  as  well  as  to 
position  this  beam  on  a  path  parallel  to  the  transmitted  beam  (bq-).  The  two 
beams  of  approximately  equal  intensity  were  then  focused  by  a  multi-element 
lens  with  a  17.4  cm  focal  length  to  form  an  interference  fringe  pattern  at  the 
sample  volume.  The  sample  volume  was  imaged  through  a  pinhole  and  then 
refocused  onto  a  fast  photomultiplier  (Type  R663)  which  views  the  light 
scattered  at  90°.  An  interference  filter  was  used  before  the  photomultiplier 
to  reject  light  outside  a  1  nm  band  centered  at  514.5  nm. 

The  measurement  volume  is  defined  by  the  intersection  of  the  laser  beam 
intensity  profile  and  the  field  of  view  provided  by  the  pinhole  of  the 
receiving  optics.  Its  geometry  is  depicted  in  Figure  5.3.  The  incident  beams 
were  at  a  7.5°  angle  with  respect  to  each  other  and  had  a  total  power  of  130 
milliwatts  with  the  iris  diaphragm  fully  open.  The  fringe  pattern  formed  from 
the  two  beams  was  easily  viewed  by  magnifying  the  measurement  volume  by  a 
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factor  of  200  with  a  microscope  lens  and  projecting  it  onto  a  screen.  The 
number  of  observed  fringes  is  consistent  with  a  calculated  4  urn  fringe  spacing 
and  12  beam  diameter  (Figure  5.4). 


A  successful  measurement  requires  the  double  coincidence  that  a  single 
particle  traverse  the  fringe  system  and  be  detected  within  the  field  of  view 
of  the  receiving  optics.  Initially,  the  measurement  volume  was  imaged  onto  a 
25  ^m  pinhole  with  a  magnification  of  unity,  but  low  counting  rates  were 
observed  with  such  a  small  field  of  view.  Also,  the  detection  of  scattered 
light  through  the  small  pinhole  was  highly  susceptible  to  random  mechanical 
vibrations.  A  400  pinhole  was  subsequently  used  and  the  signal 
detectability  was  greatly  increased,  however,  the  measurement  volume  was  too 
large.  Extraneous  signals  from  the  interaction  of  a  single  beam  with  a 
particle  arose  because  the  field  of  view  of  400  ^m  was  not  restricted  to  the 
overlap  region  of  the  two  beams.  Also,  the  larger  field  of  view  increases  the 
chance  of  detecting  multiple  particle  events.  Finally,  the  sample  volume  was 
imaged  with  a  magnification  of  =  10  onto  a  400  pinhole.  This  was 
equivalent  to  a  field  of  view  of  40  |im  and  rendered  the  signal  less 
susceptible  to  vibrations.  The  collection  optics  consisted  of  a  60  mm 
diameter  70  mm  focal  length  plano-convex  lens  located  ~77  mm  from  the 
measurement  volume.  The  overall  optics  have  a  f  number  of  ~1.3. 

5 .3  Signal  Acquisition  System 

The  data  acquisition  system  is  shown  in  Figure  5.5.  The  signal  from  the 
photomultiplier  tube  was  amplified  by  a  high  speed  buffer  amplifier  (200  MHz 
bandwidth,  xlOO  gain).  The  Brownian  motion  data  was  recorded  on  a  Data 
Precision  Flash  A/D  system  with  high  speed  memory.  For  each  event  (traversal 
of  a  particle  through  the  probe  volume),  8192  data  points  at  10  ns  interval 
were  collected.  The  record  length  of  8  ^s  of  data  was  large  compared  to  the 
relaxation  time  of  the  particle  which  is  of  the  order  of  0.1  therefore 

there  should  be  enough  data  for  statistical  analysis. 
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(a) 


(b) 


Figure  5.4. 


Magnified  View  of  the  Interference  Pattern:  (a)  With  a  Narrow 
Unfocussed  Beam  Diameter  (~0.2  um  Beam  Diameter);  (b)  With  a 
2.2  cm  Unfocused  Beam  Diameter.  The  weak  interference  fringes 
in  the  background  of  (b)  are  due  to  scattered  light  from  the 
microscopic  objective. 
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To  verify  that  the  Brownian  motion  signal  was  captured  when  the  particle 
was  located  at  the  "steep"  part  of  the  fringe  pattern,  the  signal  from  the 
buffer  amplifier  was  fed  through  a  low  pass  filter  to  provide  the  trigger 
signal  for  the  data  acquisition  system.  The  filter  frequency  was  set  at 
approximately  the  Doppler  frequency  due  to  the  convection  velocity  of  the 
particle  through  the  measurement  volume.  A  typical  trace  of  the  filtered 
signal  which  was  recorded  by  a  storage  oscilloscope  simultaneously  with  the 
high  speed  Brownian  motion  data  is  shown  in  Figure  5.6.  Because  the  trigger 
level  was  known,  the  portion  of  the  signal  captured  by  the  high  speed  data 
acquisition  system  may  be  identified  on  the  scope  trace  as  in  Figure  5.6.  If 
this  portion  of  the  signal  was  not  located  at  the  "steep"  part  of  the  fringe 
pattern  (as  shown),  the  particular  data  record  was  discarded. 

The  high  speed  data  record  corresponding  to  the  event  of  Figure  5.6  is 
shown  in  Figure  5.7.  The  top  trace  in  Figure  5.7  is  the  actual  signal  sample 
at  10  ns  Intervals.  The  bottom  trace  is  the  same  signal  with  the  points 
connected  by  a  (sin  x)/x  fitting.  Stored  data  records  similar  to  this  were 
used  for  the  statistical  analysis. 

5  .4  Results  and  Discussion 

It  was  realized  that  in  the  Initial  setup  the  signal  to  noise  ratio  was 
too  small  for  detection  of  Brownian  motion.  The  incident  beam  Intensity  of  65 
nW  focused  into  a  12  ^m  spot  (intensity  of  5  x  108  W/m4* )  was  too  low.  The 
system,  nevertheless,  was  set  up  to  verify  that  the  data  acquisition  system 
was  operational  and  to  assess  the  actual  S/N  ratio  to  determine  the 
improvements  necessary  to  obtain  a  useful  signal.  The  typical  signal  records 
collected  were  depicted  in  Figures  5.6  and  5.7.  These  records  were  processed 
digitally  for  statistical  analysis. 

To  verify  the  operation,  the  system  was  used  to  look  at  small  latex 
spheres  suspended  in  water.  Since  the  viscosity  of  water  is  '3  x  10*  times 
that  of  air,  the  signal  due  to  Brownian  motion  was  negligible  and  the  response 
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Figure  5.6. 


The  Filtered  Signal  From  the  PMT.  The  signal  was  used  to 
trigger  the  100  MHz  data  acquisition  system.  The  portion  of 
signal  recorded  Is  marked  by  the  circle  (8,192  samples  at  10 
ns  sampling  intervals). 
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Figure  5.7.  Signal  Recorded  by  the  Fast  A/D  System,  (a)  Actual  Signal 

Sample  Points  at  10  ns  Sampling  Intervals,  (b)  Sample  points 
connected  by  sin  x/x  fit. 
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of  Che  system  was  mainly  due  Co  shot  noise.  Therefore,  Che  measurement 
established  the  noise  limit  of  the  system.  The  signal  was  processed  to  obtain 
histograms  of  the  time  between  extrema  of  the  signal  as  described  in  Chapter 
3.  The  histograms  for  particles  of  several  sizes  are  shown  in  Figure  5.8, 
which  shows  that  there  was  no  difference  in  signal  response  for  different 
size  particles.  This  is  because  the  signal  was  masked  by  the  shot  noise. 

(The  peak  at  30  ns  is  probably  due  to  the  30  MHz  bandwidth  of  the  system.) 

To  improve  the  signal  to  noise  ratio  of  the  system,  two  changes  were 
made.  The  receiving  optics  were  rearranged  to  a  forward  scattering  collection 
geometry,  and  the  focusing  lens  for  the  incident  beam  was  replaced  by  one  with 
a  shorter  focal  length  (7  cm  compared  to  17.4  cm).  The  2  cm  diameter  incident 
beams  were  focused  to  a  ~4  pm  spot.  The  fringe  spacing  was  1 .5  pm  and  there 
were  3  fringes  in  the  focal  volume.  The  change  of  the  focusing  lens  improved 
the  signal  by  increasing  the  beam  intensity  and  decreasing  the  fringe 

Q 

spacing.  The  intensity  of  each  beam  at  the  measurement  volume  was  ~5  x  10 
W/m2  . 

The  above  system  was  tested  by  injecting  latex  particles  of  0.045,  0.087 
and  0.261  pm  diameters  through  the  measurement  volume.  These  particles  were 
commercially  available  as  concentrated  aqueous  suspension  in  small  vials.  The 
suspension  was  diluted  with  alcohol  and  dispersed  with  an  air  nebulizer.  The 
concentration  of  the  suspension  was  such  that  it  was  unlikely  to  have  more 
than  one  particle  in  each  nebulized  droplet.  The  droplets  were  mixed  with  a 
stream  of  drying  air  to  remove  the  liquid  content  to  ensure  that  dry  latex 
particles  were  Introduced  into  the  measurement  volume. 

The  signal  to  noise  ratio  was  calculated  from  the  measured  photocurrent 
and  the  known  gain  of  the  photomultiplier.  For  the  particles  used,  the 
photo-electron  arrival  rate  was  estimated  to  be  n  *  1010/s.  The  Brownian 

excursion  Xq,  for  these  particles  is  approximately  1  percent  of  the 
wavelength  \ .  Therefore,  for  a  30  MHz  bandwidth  system,  the  signal  to  noise 
ratio  according  to  Eq.  (5.2)  was  ~5.  Because  the  signal  processing  strategy 
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relied  on  the  frequency  modulation  (the  time  between  extrema)  to  recover  the 
necessary  information,  it  was  especially  susceptible  to  noise.  A  signal  to 
noise  ratio  of  10  is  the  minimum  acceptable  value.  Therefore  an  increase  of 
S/N  by  at  least  a  factor  of  2,  or  signal  level  a  factor  of  4  is  necessary. 

The  optical  system  had  been  designed,  approximately,  to  the  limitation  of 
commercially  available  optical  components.  The  intensity  gradient  at  the 
focal  volume  had  been  maximized  for  a  given  laser  power.  The  single  mode, 
etalon  stabilized  laser  power,  however,  may  be  easily  increased  by  using  a 
newer  model  laser.  Alternatively,  the  instrument  may  be  used  to  look  at 
particles  in  the  >  0.2  |j,m  range  (the  S/N  is  proportional  to  db).  At  a  higher 
beam  intensity,  however,  other  physical  processes  begin  to  play  an  important 
role  in  modifying  the  Brownian  motion  behavior.  These  processes  and  the 
limitations  imposed  by  them  are  discussed  in  the  next  chapter. 


6.  DETECTION  LIMIT  OF  BROWNIAN  MOTION 


The  scattered  light  from  a  particle  may  be  increased  by  increasing  the 
intensity  of  the  probe  laser  beam.  There  is  a  limit,  however,  on  the  maximum 
beam  intensity  because  it  will  disturb  the  dynamics  of  the  particle  Brownian 
motion.  There  are  two  effects:  (a)  the  particle  may  be  heated  so  that  the 
thermal  environment  changes,  and  there  is  a  substantial  thermophoretic 
effect.  (We  are  not  considering  the  severely  heated  case  in  which  the 
particle  may  vaporize  or  Ignite.)  (b)  At  high  beam  intensity,  the  photon 
pressure  may  change  the  dynamics  of  the  particle  motion.  The  thermophoretic 
effect  depends  on  the  particle  absorptivity  and  the  thermal  properties  of  the 
particle  and  the  fluid.  The  photon  pressure  effect  is  present  Independent  of 
the  particle  absorptivity.  These  effects  determine  the  limit  of  detection  of 
the  Brownian  motion.  In  the  following,  the  above  effects  will  be  examined. 

The  particles  are  assumed  spherical  and  the  diameters  of  interest  are  in  the 
range  of  0.02  ^.m  to  0.2  ^m.  A  diameter  of  0.1  is  used  as  the  nominal  value 
to  estimate  the  effects  quantitatively.  The  laser  beam  power  P  is  assumed  to 
be  1W  at  0.5  |io  wavelength  and  is  focused  into  a  spot  of  diameter  b  -10  o» . 

The  beam  Intensity  is  ~1010  W/ta"  .  These  values  would  correspond  to  a  1  cm 
diameter  beam  focused  to  a  spot  by  a  10  cm  focal  length  convex  lens.  The 
thermal  environment  of  the  particle  is  atmospheric  air  at  300  K.  The 
particles  are  assumed  to  be  convected  through  the  probe  volume  at  a  velocity  U 
-  1 0  cm/ s . 

Since  the  mean  free  paths  A  of  the  air  molecules  at  atmosphere  is 
-0.06  um,  particles  in  the  range  of  0.02  to  0.2  am  are  in  the  transition 
flow  regime  (Knudsen  number  kn  ■  A/d  ~3  to  0.3).  The  theory  of  transition 
flow  has  not  been  well  established.  The  present  calculations  were  based  on 


6.1  Thermal  Characteristics  of  the  Particles 


In  the  flow  regime  considered  here,  heat  transfer  from  the  submicron 
particles  by  forced  convection  and  natural  convection  are  negligible  due  to 
the  extremely  small  Reynolds  number  and  Grashof  number.  The  major  heat 
removal  mechanism  is  through  conduction.  Data  on  the  heat  transfer  from  a 
sphere  in  a  rarefied  gas  from  Ref.  [6.2]  may  be  fitted  as: 


U  ,  ( 

—  -  1  -exp( 


0.215  0.0275  _  2.815  x  10 

2  3 

k  k  k  J 

n  n  n 


1.346  x  10 


2.295  x  10 


n  (6.1) 


where  the  Nusselt  number  (Nu)  is  evaluated  at  the  free  stream  properties  and 
Nuc  is  the  Nusselt  number  at  the  continuum  limit  (Nuc  -  2).  With  a  0.1  pm 
diameter  particle  in  atmospheric  condition  *  Nu  ~0.5.  Values  for  Nu  as  a 
function  of  kn  are  shown  in  Figure  (6.1). 

As  the  particle  enters  the  beam  volume,  it  is  heated  up  by  the  laser  beam 
because  of  the  finite  absorptivity.  The  heat  up  time  may  be  estimated  by 
considering  the  thermal  diffusion  time  through  the  particle  as  well  as  the 
boundary  layer.  For  an  order  of  magnitude  estimate,  we  used  the  properties  of 
polyethylene  particles  with  a  thermal  diffusivity  ap  ~8  x  10“  m  / s.  The 
temperature  distribution  inside  the  particle  will  become  steady  state  in 
tp  ~  d2/<ip.  For  particles  in  the  0.02  to  0.2  ^m  range,  tp  ~0.005  to 
0.5  ^18.  The  thermal  diffusion  time  through  the  boundary  layer  may  be  obtained 
from  the  energy  balance  on  the  particle 


me  77  ■  la  - N  xd2  k(T-T  )/d 

p  dt  4  u  <*> 


(6.2) 


where  m,  cp  are  the  mass  and  specific  heat  of  the  particle,  I  Is  the  laser 
beam  Intensity,  a  the  absorptivity  of  the  particle,  and  k  the  thermal 
conductivity  of  the  air.  The  thermal  diffusion  time  from  Eq.  (6.2)  is 
therefore : 
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The  steady  state  temperature  rise  is  plotted  against  the  product  la  in 

Figure  6.2  for  particles  in  the  0.02  to  0.2  pm  range.  To  avoid  a  substantial 

7  2 

perturbation  of  the  thermal  environment,  la  should  be  less  than  10  W/m  .  For 


10  2 

the  1W  beam  focused  to  10  pm  square  beam  spot  (1  ~10  W/m  ),  the  particle 


absorptivity  should  be  10“3  or  less.  Thus  the  method  cannot  be  used  for 
particles  such  as  soot  particles  which  have  a  high  absorptivity. 


6 .2  Thermophoresis  of  the  Particles 


When  a  particle  with  non-zero  absorptivity  is  illuminated  by  a  laser 
beam,  thermophoresis  occurs  because  a  temperature  gradient  is  set  up  so  that 
the  illuminated  side  of  the  particle  has  a  higher  temperature  than  the  shadow 
side.  The  net  result  is  that  the  particle  acquires  a  drift  velocity  in  the 
laser  beam  direction.  This  problem  may  be  studied  in  a  similar  manner  to  the 
classical  analysis [6. 3]  of  thermophoresis  of  small  particles  in  thermal 
equilibrium  with  a  gas  in  which  a  temperature  gradient  is  present.  In  the 
present  case,  however,  the  boundary  condition  on  the  particle  surface  would 
have  to  be  modified  to  include  the  heat  flux  due  to  the  absorption  of  the 
laser  light,  and  the  temperature  field  far  from  the  particle  would  be 
uniform.  Such  a  case  has  not  been  discussed  in  the  literature.  To  develop 
the  solution  for  this  case  would  be  beyond  the  scope  of  this  study.  We, 
therefore,  use  a  simple  one  dimensional  energy  balance  estimate  of  the 
temperature  gradient  across  the  particle,  and  then  estimate  the  thermophoretic 
effect  using  this  temperature  gradient  and  the  semi-empirically  based  formula 
suggested  by  Ref.  [6.4]. 


In  the  one  dimensional  model  of  this  "spherical"  particle,  the  particle 
is  modelled  as  a  cylinder  of  diameter  d  and  length  d  as  shown  in  Figure  6.3. 
The  side  wall  of  the  cylinder  is  assumed  to  be  adiabatic.  Each  end  of  the 
cylinder  is  fitted  with  a  heat  exchanger  of  area  nd2/2,  representing  the 
surface  area  of  one  half  the  original  sphere.  The  heat  transfer 
characteristics  of  these  heat  exchangers  are  specified  through  the  Nusselt 
number  of  a  sphere  in  a  rarefied  flow.  The  temperature  profile  is  shown  in 
Figure  6.3.  The  power  input  from  the  laser  beam  is  la  ud  /4.  In  steady 
state,  therefore, 
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The  average  temperature  Tave  =  (Tj  +  T2)/2  and  the  temperature  difference  AT 
=  Tj  -  T2  may  be  calculated  from  Eq.  (6.4)  and  Eq.  (6.5)  as 


T 

ave 


-  T 

00 


la  d 
4  Nu  k 


(6.6) 


la  d  v  ave 

T  2^kp  +  Nu  k)  [l  +  kp/(k  Nu) j  (6*7) 

The  temperature  difference  AT  la  plotted  In  Figure  6.4.  It  Is  a  significant 
fraction  of  the  average  temperature  rise  (Tave  -  T„). 

In  the  near  continuum  limit,  the  thermophore  tic  force  F  for  a  particle  in 
equilibrium  with  a  gas  In  which  a  temperature  gradient  Is  present  was 
calculated  by  Brock[6.3] 


2 

6wU  d  C  [k/k  +  2  C  A/d 

p  _ _ 8  'VT'|  P  t _ 

oT.  '  x  [1  +  6  C  A/d  1  +  2  k/k  +~UC  A7d^ 
u  m  p  t 

where  C8  (-  1.17)  Is  the  thermal  slip  coefficient,  Ct  (-2.18)  Is  the 
temperature  jump  coefficient,  and  Cra  (-1.14)  Is  the  momentum  exchange 
coefficient.  Talbot[6.4]  argued  that  this  formula  is  a  useful  means  of 
interpolation  for  all  values  of  A/d,  since  It  reduces  to  within  of  the 
correct  colllslonless  limit  (for  perfect  accommodation)  when  A/d  *  ». 
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Under  the  action  of  the  thermophore tic  force,  the  particle  would  attend  a 
terminal  velocity  in  a  time  scale  of  the  order  of  the  Brownian  relaxation  time 
(relax  =  pp  d2/18  p).  Using  the  Milliken  drag  formula, [6 .4 ]  the  drift 
velocity  is 

2vC  (k/k  +  2C  A/d')  fl  +  ~  f A  +  B  e~Cd/2A)l  ?T 

UT  "  (1  +  60  A/d)  (1  +  2  k/k  +  4  C  A/d)  ri7_^  (6,9) 

v  m  '  K  p  t  0 

where 

A  *  1 .20 
B  -  0.41 
C  =  0.88 


Although  the  assumptions  leading  to  Eq.  (6.9)  do  not  match  the  present 
application,  (most  importantly,  the  temperature  distribution  around  the 
particle  would  be  quite  different),  a  useful  estimate  of  the  drift  velocity 
may  be  obtained  by  replacing  the  temperature  gradient  in  Eq.  (6.9)  with  AT/d, 
with  AT  given  by  Eq.  (6.7).  If  the  temperature  rise  AT  is  not  very  large, 

To,  may  be  used  for  T0 .  The  the rraopho retie  drift  velocity  of  the  particle 
and  the  mean  Brownian  velocity  [/(8  kT/nrnp)],  are  plotted  against  the 
particle  size  in  Figure  6.5.  The  thermophoretic  velocity  uq-  is  not  very 
sensitive  to  the  particle  diameter  because  AT  is  proportional  to  d  so  that 
(VT)X  is  Independent  of  d.  At  a  modest  absorbed  power  density  level  of  107 
W/m2,  at  which  the  temperature  rise  of  the  particle  is  ~20°C  according  to 
Figure  6.2,  the  drift  velocity  is  significantly  higher  than  the  mean  Brownian 
velocity.  Therefore  thermophoretic  effects  could  significantly  alter  the 
statistics  of  the  Brownian  velocity. 


6-9 


6.3  Photon  Pressure 
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At  high  laser  beam  intensity,  there  is  a  significant  pressure  force  on 
the  particle  due  to  the  momentum  exchange  between  the  photons  and  the 
particle.  Under  the  motion  of  this  force,  the  particle  will  attend  a  terminal 
velocity  in  a  time  scale  of  the  order  of  the  Brownian  relaxation  time 

a 

(x relax  ~  pp  d  /18  p).  If  the  laser  beam  is  uniform  and 

uni -directional,  this  terminal  velocity  would  not  affect  the  Brownian  motion 
statistics.  The  actual  beam  probe  volume,  however,  is  obtained  by  focusing 
the  laser  beam  to  a  diffraction  limited  spot.  As  a  result,  both  the  intensity 
and  direction  of  the  laser  light  is  highly  dependent  on  the  position  of  the 
particle,  and  the  photon  pressure  effect  is  coupled  into  the  Brownian  motion. 
To  estimate  the  effect,  we  shall  calculate  the  terminal  velocity  of  the 
particle  due  to  the  photon  pressure  and  compare  it  to  the  mean  Brownian  motion 
velocity. 


Consider  a  particle  which  is  illuminated  by  a  polarized  laser  beam  in  the 
z  direction  as  shown  in  Figure  6.6.  The  incident  electric  field  is  in  the  x 
direction.  According  to  Mle  theory, [6.5]  the  electric  field  intensity  of  the 
scattered  light  is  given  by 


v, 


/i^e)  o 


0  /i2(ell  L 


(6.10) 
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where  the  scattering  amplitudes  ij  and  i2  are  in  terms  of  the  size  parameter 
x(-rtd/\)  and  the  refractive  index  m  of  the  particle.  If  I0  is  the  incident 
beam  intensity  ( photons/area/tirae) ,  then  the  force  in  the  z  direction  due  to 
the  momentum  transfer  of  the  scattered  photon  is 


lQ  h  >  *  2  2 

F  »  — r  —  j  d(j>  /  de  sin  9  (i-cos  9)  Ci  (9)  sin  $  +  i  ( © )  cos  $1 
2  kZ  X  0  0  1  2  (6.11) 
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where  h  is  the  Planck's  constant  and  the  $  dependance  is  illustrated  in 
Figure  6.6.  Since  k  *  2 n/X,  and  i^  ,  i2  are  only  functions  of  cos  9,  (see 
Ref.  [6.5]),  the  force  due  to  photon  pressure  is 

1 

-1 

cos  0)  [i,(co8  9)  +  i_(cos  0)J  (1-cose) 

1  2  (6.12) 

where  I  is  now  the  beam  intensity  in  W/m2 .  To  obtain  representative  values  of 
Fz,  the  integral  in  Eq.  (6.12)  was  evaluated  numerically.  For  a  0.25  jim 
particle  with  refractive  index  i  «  1.5  +  0  1,  which  is  representative  of  the 
polystyrene  test  particles,  the  values  for  Fz  are  shown  in  Figure  6.7.  Also 
shown  are  the  corresponding  drift  velocities  uj)»  given  by  the  Milliken 
formula: 


d(cos  9)  [i^(cos  9)  +  ^(cos  9)]  (l-cos9) 
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(6.13) 


The  results  show  that  for  non-absorbing  particles  of  0.25(im  diameter,  the 
photon-pressure-produced  drift  velocity  is  less  than  the  Brownian  velocity  for 
beam  intensity  up  to  10AU  W/m2.  Since  the  scattering  cross  section  decreases 
rapidly  with  particle  size,  photon  pressure  would  not  significantly  affect  the 
Brownian  motion  for  beam  intensities  of  <10io  W/m2. 
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6.4  References  for  Chapter  6 


7.  SUMMARY  AND  CONCLUSIONS 


A  cheoredcal  and  experimental  study  on  the  feasibility  of  determining 
the  size  of  subalcron  particles  by  observing  their  Brovmlan  motion 
characteristics  has  been  carried  out.  The  motion  was  detected  by  using  an 
interferometric  system  very  similar  in  construction  to  a  conventional  laser 
doppler  veloclmeter  system.  The  mean  excursion  associated  with  the  Brownian 
motion,  however,  is  much  smaller  than  the  fringe  spacing  formed  by  the 
intersecting  laser  beams  at  the  measurement  volume  so  that  it  is  not  possible 
to  extract  information  from  the  signal  using  conventional  laser  doppler 
veloclmeter  signal  processing.  Nevertheless  it  is  possible  to  determine  the 
damping  behavior  of  the  particle  in  Brownian  motion  by  analysing  the  frequency 
content  of  the  signal.  This  damping  is  characterized  by  a  relaxation  time 
which  is  a  function  of  the  temperature  of  the  fluid  in  which  the  particle  is 
Immersed,  the  viscosity  of  the  fluid,  and  the  inertia  of  the  particle.  If  the 
fluid  temperature  and  viscosity  are  known,  the  Inertia,  and  thus  the  size  of 
the  particle,  may  be  determined  from  the  relaxation  time  which  is  obtained 
from  the  Brownian  motion  detector. 

A  theoretical  analysis  of  the  signal  processing  strategy  was  conducted 
using  a  Monte  Carlo  simulation  of  the  particle  motion.  In  the  simulation  it 
was  possible  to  determine  the  relaxation  time  of  the  particle  from  the  mean 
value  of  the  time  between  the  extrema  of  the  signal,  which  represented  the 
zeros  of  the  Brownian  velocity  fluctuation.  Preliminary  analysis  of  the 
design  of  the  experiment  indicated  that  Che  signal  to  noise  of  the  Brownian 
detection  system  was  marginal.  The  performance  of  the  system  depended  very 
much  on  the  actual  Implementation  and  had  to  be  evaluated  experimentally. 

A  proof-of-concept  experiment  was  conducted  using  submicron  latex  spheres 
of  known  diameters  as  test  particles.  The  optics  had  been  designed  to 
approximately  the  limit  of  commercially  available  optical  components.  The 
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laser  power,  however,  was  not  adequate  to  achieve  the  required  S/N  for 
retrieving  the  statistical  Information  on  the  particle  Brownian  motion.  (The 
argon  ion  laser,  which  nominally  puts  out  2.5  W  in  a  single  line,  was  found  to 
lose  75  percent  of  its  power  when  operated  with  a  mode  selection  etalon  which 
was  necessary  for  the  present  experiment.)  The  recorded  signal  was  found  to 
be  dominated  by  shot  noise.  It  was  estimated  that  the  laser  power  density 
needed  to  be  increased  by  at  least  a  factor  of  4  to  achieve  adequate  signal  to 
noise  ratio. 

While  the  required  laser  power  density  may  be  easily  obtained  by  using  a 
newer  model  laser,  other  physical  processes  begin  to  influence  the  Brownian 
motion  characteristics  at  high  laser  beam  intensity.  An  analysis  of  the 
effects  of  thermophoresis  and  photon  pressure  was  carried  out.  The  effect  of 
thermophoresis  due  to  the  uneven  heating  of  the  particle  by  the  laser  beam  was 
found  to  be  a  major  limitation  so  that  the  application  of  the  method  to 
particles  with  finite  absorptivity  (greater  than  ~10-4)  may  lead  to 
erroneous  results.  The  effect  of  photon  pressure  was  found  to  be  significant 
for  beam  intensity  greater  than  I010  W/m2 . 

For  an  experiment  to  demonstrate  the  principle,  particles  of  low 
absorptivity  (such  as  polystyrene)  may  be  used.  Our  present  results  suggest 
that  for  particles  In  the  0.1  u®  range  there  would  be  adequate  S/N  (~20)  if 
a  stable  single  mode  laser  with  power  of  2.5  W  is  used. 

Although  the  application  to  soot  particles  originally  envisioned  will 
encounter  the  limitations  on  laser  beam  Intensity  described  above  (see 
Subsection  6.1),  other  applications  may  exist,  such  as  refractory  oxides  in 
relatively  clean  combustion  systems. 
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8.  ADMINISTRATIVE 


8.1  Personnel 

The  Principal  Investigators  for  this  work  through  almost  the  entire 
contract  period  were  Dr.  Alan  Stanton,  a  Principal  Research  Scientist  at 
Aerodyne  Research,  Inc.  and  Dr.  Wai  Cheng,  Principal  Research  Scientist 
(part-time)  at  Aerodyne  and  Associate  Professor  of  Mechanical  Engineering, 

MIT.  Dr.  Stanton  concentrated  his  research  efforts  for  this  program  on  the 
development  of  optical  and  data  acquisition  systems  for  the  experimental 
program,  while  Dr.  Cheng  concentrated  on  the  measurement  theory,  through 
development  of  the  Monte  Carlo  Brownian  motion  simulation  model.  In  addition 
to  his  work  at  Aerodyne  on  this  project,  Dr.  Cheng  supervised  the  work  of  two 
graduate  students  in  the  Department  of  Mechanical  Engineering,  MIT,  Sarah 
Kostic  and  Eugenie  Halnsworth.  This  work  was  in  support  of  the  Brownian 
motion  simulation  model,  funded  by  a  subcontract  to  MIT. 

The  laboratory  implementation  of  this  program  was  mostly  performed  by 
Mr.  Keith  McCurdy,  who  joined  the  staff  of  Aerodyne  in  the  second  year  as  a 
Research  Scientist  after  receiving  his  Masters  degree  in  Physical  Chemistry 
from  Rice  University. 

In  January  of  1986  Dr.  Stanton  left  Aerodyne,  and  his  duties  as  Principal 
Investigator  were  assumed  by  Dr.  Joda  Wormhoudt.  Dr.  Cheng  continued  as 
Principal  Investigator  at  MIT  and  supervised  the  final  experimental  work. 

8 .2  Interactions 

Three  presentations  were  made  at  AFOSR  Contractors'  Meetings  on  this  work 
during  the  first  year  of  the  contract,  one  during  the  second  year,  and  one 
during  the  third.  These  presentations  were: 


First  Year 


A.C.  Stanton  and  W.K.  Cheng,  "Techniques  for  Subraicron  Particle  Sizing," 
1983  AFOSR  Contractors  Meeting  on  Air  Breathing  Combustion  Dynamics 
Research,  Scottsdale,  AZ,  September  19-22,  1983. 

W.K.  Cheng  and  A.C.  Stanton,  "Single  Particle  Sizing  by  Measurement  of 
Brownian  Motion,"  1984  AFOSR  Research  Meeting  on  Diagnostics  of  Reacting 
Flows,  Yale  University,  March  21-22,  1984. 

A.C.  Stanton  and  W.K.  Cheng,  "Single  Particle  Sizing  by  Measurement  of 
Brownian  Motion,"  1984  AFOSR/ONR  Contractors  Meeting  on  Air  Breathing 
Combustion  Research,  Pittsurgh,  PA  June  10-21,  1984. 

Second  Year 

W.K.  Cheng,  A.C.  Stanton,  and  K.  McCurdy,  "Single  Particle  Sizing  by 
Measurement  of  Brownian  Motion,"  1985  AFOSR/ONR  Contractors  Meeting  on 
Particle  Emission  Technology,  Monterey,  CA,  April  16-18,  1985. 

Third  Year 

W.K.  Cheng,  K.  McCurdy,  and  J.  Wormhoudt ,  "Study  of  Submicron  Particle 
Size  Distributions  by  Laser  Doppler  Measurement  of  Brownian  Motion,"  1986 
AFOSR  Contractors  Meeting  on  Diagnostics  of  Reacting  Flows,  Stanford,  CA, 
June  16-17,  1986. 


8.3  Publications  and  Patents 


No  publications  or  patents  resulted  from  the  work  reported. 
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